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Abstract 

A simple deterministic model is proposed to represent the ba- 
sic aspects concerning the effects of different antiretroviral treatment 
schedulles on HIV incidence and prevalence of affected populations. 
The model mimics current treatment guidelines applied in Brazil. 
However, the model does not intend to fit the data with any acceptable 
degree of accuracy since uncertainties on the values of the parameters 
and on the precise effect of the treatment put some limits on the prac- 
tical implications of our model from which only orders of magnitude 
and some qualitative effects can be deduced. So, this paper intends to 
provide a conceptual and mechanistic understanding of the possible 
long term effects of treatment on the dynamics of HIV transmission. 

According to the model, the effect of the treatment depends on the 
level of sexual activity of the subpopulations considered, being more 
pronounced on the subpopulations with the highest sexual activity 
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levels. Also, inefficient treatment can be prejudicial depending on 
the level of sexual activity and on the capacity to provide adequate 
treatment coverages to the population affected. 

1 Introduction 

The natural history of the Human Immunodeficiency Virus (HIV) infection 
and transmission is well understood today. It is now accepted that the viral 
concentration in the circulating blood (viraemia) and other organic fluids 
determines the probability of transmission ([1. It has been demonstrated by 
several authors that after infection the viraemia rapidly increases in the first 
week to a high level, staying in that high levels for some weeks, dropping 
afterwards to very low levels and then increasing slowly for several years. 
After a period of 10 - 15 years, the increasing level of viraemia reaches a level 
in which clinical manifestations begin, coinciding with the dropping in the 
CD4 lymphocytes (the target cell of HIV) counting and the development of 
the conditions that define the state of full-blown AIDS. 

So, it should be expected a large variation in the probability of transmis- 
sion of HIV along the natural course of the infection [2]. The relative con- 
tribution of each of these distinct phases of viraemia seen along the natural 
history of the infection to HIV transmission have only recently been demon- 
strated jSj- This relative contribution may have important consequences to 
the epidemiological pattern of HIV transmission, as well as to the assess- 
ment of the impact of HIV treatment on the epidemiology of HIV/AIDS. In 
addition, it should be expected that differences in the probability of transmis- 
sion in each of the distinct phases of the infection may also have important 
consequences on the evolution of HIV virulence [3] . 

The use of the combined antiretrovirus therapies (ARVTs), particularly 
those known as HAART (highly active antiretrovirus therapies) has demon- 
strated excellent results in all the countries where people infected by HIV 
has access to the treatment 0. The results reported significant reduction in 
deaths rates caused by opportunistic infections. For example, in the U.S.A. 
there has been a reduction of 61% between 1995 and 1997 in the mortality 
due to AIDS. In Brazil, whose experience has been acclaimed worldwide as 
one of the most successful attempts to control AIDS, the government de- 
cided to provide ARVT to all HIV seropositive individuals who fulfilled the 
treatment criteria. From 1996 onwards HAART has been used as a standard 
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treatment, with 120,000 patients being treated in 2002. 

The treatment against HIV evolved rather rapidly. Just 4 years after the 
identification of HIV as the causative agent of the new syndrome, AIDS, zi- 
dovudine (AZT) the first drug for its treatment was licenced by the FDA. In 
the next eight years three other nucleoside analogs (the same class as zidovu- 
dine) were introduced. Concurrently, a better understanding of the dynam- 
ics of HIV replication and drug resistance mechanisms caused a shift from 
single- to combination-drug therapy. From 1995 to 1998 eight new antiretro- 
viral agents were approved, including protease inhibitors and non-nucleosides 
transcriptase inhibitors . These latter drugs also provided clinicians with a 
highly effective antiretroviral therapy, known as HAART that reduces viral 
load by a factor of 10 3 . However, these new drugs cause several important 
side effects that limit their use. Because of this, current clinical guidelines 
propose a delaying in starting HIV treatment to a point in which either 
the viral load is high (above 30000 copies per ml) or the CD4 cell count is 
low. However, with the progression of the infection, even in the persistence 
of treatment, HIV viral load tends to increase again leading to a need to 
change the drugs in use or to progression to AIDS [U]. 

The discussion on the implications of antiretroviral antiretroviral treat- 
ment has been restricted to clinical and virological aspects. However, con- 
sidering that different antiretroviral treatment strategies have distinct effects 
on viral load, one crucial aspect to be considered is the epidemiological con- 
sequences of a given antiretroviral treatment strategy on the incidence of 
new HIV infections. On one hand, effective antiretroviral treatment reduces 
viral load. On the other hand, it prolongs the asymptomatic phase, proba- 
bly the most important for transmission [3| in the absence of antiretroviral 
treatment. Both effects are intrinsically related to transmission, although in 
opposite directions. Therefore, the critical question related to the choice of 
the best antiretroviral strategy should take into account the epidemiological 
consequences of different antiretroviral treatments. In spite of some attempts 
to understand this crucial aspect, the best antiretroviral treatment strategy 
is still to be defined [7|,[E], 0. 

In this paper we propose a schematic mathematical model to analyze the 
impact of current antiretroviral therapy on the incidence and prevalence of 
HIV infection. In section 2, we point out some historical aspects of Brazilian 
public health strategies related to HIV/AIDS treatment and its impact on the 
epidemics in Brazil. Our model was constructed incorporating the treatment 
concepts of the Brazilian program. However, due to uncertainties in param- 
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eters, practical conclusions should be taken with caution. In section 3, we 
derive an integral equation which allows the calculation of the impact on HIV 
incidence of an antiretroviral treatment program. This equation depends on 
a number of functions and parameters whose forms and values are described 
in section 4. In section 5, we present numerical results of several hypothetical 
scenarios and in section 6 we discuss the model's limitations and draw some 
tentative conclusions regarding possible epidemiological implications of our 
results. 

2 Motivation: the Brazilian experience with 
HAART 

Since 1991 the Brazilian government decided to provide ARVT to all HIV 
seropositive individuals who fulfilled the treatment criteria ^H]- From 1996 
onwards HAART has been used as a standard treatment, with 120,000 pa- 
tients being treated in 2002. Its impact was immediately noted by health 
authorities, in particular its dramatic reduction on AIDS mortality, and this 
public health experience has been acclaimed worldwide as one of the most 
successful attempts to control AIDS. 

Figure 1 shows the reduction in the mortality by AIDS in the period 
between 1996 and 2001. The total number of averted deaths summed up to 
90,000 patients in this period. 

Figure 1 

In addition, it was observed an increase in the survival period of patients 
with AIDS, from 5 to 58 months, a 12 fold increase, after the introduction 
of HAART treatment. Also, there was a six- fold reduction in the number of 
AIDS hospitalizations and a reduction of 54% in the cost of treatment. 

With respect to the impact of HAART on the incidence and prevalence 
of HIV infection, only rough estimates are available. For instance, in 1992 
the World Bank projected the number of expected HIV infected individuals 
in Brazil by 2002 as 1.2 million. However, the current estimates for 2002 
is around 600 thousands. It is difficult to attribute this reduction solely 
to HAART, since other preventive practices have been highly stimulated 
by the Brazilian government. For instance, the observed reduction in HIV 
infection prevalence among sex workers from 18% in 1996 to 6.1% in 2000, and 
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among homosexuals from 10.8% in 1999 to 4.7% in 2001, has been attributed 
by health authorities to the association of condom distribution and other 
preventive measures adopted. 

The most important indication of the possible impact of HAART on HIV 
transmission in Brazil is the marked reduction in incidence rates of AIDS, 
as shown in figure 2. It is noticeable that after a historical increase in the 
incidences from 8 cases per 100,000 people per year in 1991 up to the peak 
of 18 cases per 100,000 people per year attained in 1998, the number of new 
cases dropped to 5 cases per 100,000 people per year in 2003. If we consider 
that the universal treatment started in 1996, we may conclude that after a 
delay of two years, the effect of HAART on the incidence of AIDS could be 
observed. 

Figure 2 

Notwithstanding, the observed reduction in HIV incidence and prevalence 
cannot be entirely explained by the effects due to HAART. Additionally, 
behavioral changes attributable to the preventive campaigns carried out in 
Brazil simultaneously to the beginning of HAART may also have had a pos- 
itive impact on HIV transmission. For instance, there have been a marked 
increase in regular condom use verified between 1999 and 2000 (42% to 64%). 

The Brazilian experience motivated us to model the possible role of ARVT 
on the incidence and prevalence of HIV. However, lack of data and uncertain- 
ties on the precise effect of the treatment put some limits on the practical 
implications of our model from which only orders of magnitude and some 
qualitative effects can be deduced. 

3 The model: formalism 
3.1 Sexually transmitted HIV 

As in previous papers, we consider a community of N individuals in steady 
state with respect to time j^j, divided into classes according to the contact 
pattern and transmission intensity of HIV JT] . For simplicity, in this paper 
we consider that the interaction between individuals from different classes 
is so rare that can be neglected. The classes will be described in section 5. 
As the interclass interactions are neglected, we describe a general formalism 
which applies to any and all classes. 
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Let N(a)da be the number of individuals with age between a and a + da. 
Among those, let X(a)da be the number of susceptibles with age between a 
and a + da. The number of unprotected sexual contacts per unit time those 
individuals with age between a and a + da make with all other individuals 
with ages between a' and a' + da' is: 

(3 (a,a')da'X(a)da (1) 

Let Yx(a! ,r)da'dr be the number of individuals with age between a' and 
a' + da' infected when aged between r and r + dr. The number of unprotected 
sexual contacts per unit time those individuals with age between a and a + da 
make with all other individuals with ages between a' and a' + da' is 

n i i\ i ,Yi(a',T)da'dT 

(3 {a,a)da — — X(a)da (2) 

JM (a')da' 

Let g(a' — r) be the probability of a susceptible individual to get infected 
when making an unprotected sexual contact with an individual Yi(a', r)da'dr. 
This depends on viral load which in turn depends on the time since infection 
(a' — r), and its form is given by equation (jHBj) Therefore, the number of new 
infections per unit time, due to individuals Yi(a', r)da'dr, is given by 

[3 (a, a')g(a' — r) j^ryr da'drX(a)da (3) 

Let Y 2 (a' ,T,l)da'drdl be the number of individuals with age between a' 
and a'+da' infected when aged between r and r+dr and treated continuously 
after /, that is, the treatment is initiated between / and / + dl. Therefore, 
the fraction of the unprotected sexual contacts given by equation ((TJ), with 
individuals with age a' and a' + da' infected when aged between r and r + dr 
and treated between I and I + dl is 

n ( i\ i ,Y 2 {a',T,l)da'drdl 

(3 (a } a)da — X(a)da 4 

JM {a'jda' 

Note that we assumed that the sexual behavior of susceptible towards 
treated and untreated infected individuals is the same. 

Let g\ (a' — 1,1 — t) be the probability of a susceptible individual to get in- 
fected when making an unprotected sexual contact with a treated individual 
Y 2 (a', r, l)da'drdl. This depends on viral load which in turn depends on the 
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time since infection (a' — r) and the time since the start of the treatment, and 
its form is given by equation (J38|) . Therefore, the number of new infections 
per unit time, due to treated individuals Y 2 (a' ,r,l)da'drdl, is given by 

/3 (a, d)gi(a' — 1,1 — r) k ~ da'dTdlX(a)da (5) 

After integrating equation Q from to a' with respect to r, and from 
to infinite with respect to a', and integrating (j5J) from r to a' with respect to 
I, from to a' with respect to r, and from to infinity, with respect to a', 
and summing the two contributions, we get the so-called per capita force of 
infection, A (a), defined as 

A(a) = f °° f a ' (3 (a,a')g(a' -r)^^drda' 

(6) 

+ /o°° lo if A>(a, -U- T)^0rdldTda> 

As we showed in a previous paper [T2], the contact function (3 (a, a') must 
satisfy a symmetry relation. To see this, let 

P Q (a,a')da'N(a)da (7) 

be the number of unprotected sexual contacts individuals with age between 
a and a + da make with all individuals aged between a' and a' + da'. This 
number should be equal to the number of unprotected sexual contacts that 
individuals aged between a' and a' + da' make with all individuals with age 
between a and a + da, 

(3 (a',a)daN(a')da' (8) 

that is, 



N(a>) N{a) [ ' 

This can be satisfied if /3 (a, a') is of the form 

P Q (a,a')=f(a,a')^- (10) 

where A^ = J °° N(a)da is the total population and f(a,a') is a symmetric 
function of a and a', describing the per capita rate of unprotected sexual 
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contacts(see equations I2H1 and EOJ) . Note that, due to the division by popula- 
tion size an increase in the frequency of age groups entirely uninterested for 
a given person will decrease his/hers sexual activity with those specific age 
groups. However, this is reflected automatically in the form of f(a, a'). 
Substituting (fTUj) in (0) we get 



A(a) 



I roo ra 

N Jo Jo 

I roo ra' ra 

N Jo Jo Jt 
Y 2 (a', r, l)dldrda 



+ 



f(a, a')g(a' — r)F 1 (a / , r)drda' 
f(a,a')gx(a' -1,1-t) 



Now, the equation for X(a) is 

dX(a) 
da 



X(a)X(a) — fj,X(a) 



(11) 
(12) 



where \x is the natural mortality rate of humans. 
Equation (fT2j) can be integrated, resulting in 



X(a) = X(0)exp 



X(s)ds — fia 



(13) 



Let us define hi (a, r) as a function describing the removal of individuals 
from the first infective condition by natural mortality and additional mortal- 
ity due to progression to AIDS (see equation l33j). and h 2 (a, r) describing the 
removal by antiretroviral treatment (see equation 13 7|). Then we can write: 



Yi(a, r) = Yi(r, r)hi (a, r) h 2 (a, r) 



Now we have: 



Yi(a,a) = X(a)X(a) 
Substituting equation (JT3~j) in (JT3J), we have: 



rx(a,a) = X(0)A(a)exp 

Substituting (jTfij) in (fTlj) we get 

Ki(a,r) = X(0)A(r)exp 

hi (a, r) h 2 (a, r) 



— / X(s)ds — fj,a 
Jo 



(14) 



(15) 



(16) 



X(s)ds — jir 



(17) 
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Substituting equation (|T7jl in (jTTjl we get: 
X(0) r°° r a ' 



N Jo 

Xi^ei-Io^-^h, (a',r) h 2 (a',r) c/rda' 
1 

+ 7v 



JO 



/(a, a') 

#i(a' - /, / - r)Y 2 (a', r, l)dldrda' (18) 
Assuming a antiretroviral treatment rate i/(r, a), we have 

F 2 (/,r,/) = F 1 (/,r)z/(r,0 (19) 

and 

y 2 (o,r,0 = Y 2 {l,r,l)h 3 {a,r,l) (20) 

where /13 (a, r, /) is a function describing the removal of individuals from the 
treated condition by mortality (see equation EH below) . Then, substituting 
equations (fP7j) in (fTUj) and the resulting equation in equation (|2*UJl . we get 



/ \(s)ds — /xr 

JO 



K 2 (a,r,f) = X(0)A(r)exp 

h x {l,r) h 2 {l,r)h 3 {a,r,l)u{r,l) (21) 

so that we have 

A(a) = f^ (T )e-i;^ 
AT jo jo 

/(a, a')g(a' — t)H\ (a 1 , r) /i 2 (a', r) drda! 

X(0) r°° r a ' r a ' 
+ nI I I v ^ l ) h ^r)h 2 {l,T)h^{a , ,r,l) 

f(a,a')gi(a — 1,1 — r)dldrda' (22) 

Equation (|22|) always has A (a) = as a solution. Depending on the 
functions f(a, a')g(a' — r), f(a, a')g\{a! — 1,1 — r) and on the parameters of 
/(a, a')g(a' — r), /(a, a')gi(a' — 1,1 — r), /ij (a', r), /i 2 ( a ', r ), ^3 ( a ', T , an d 
z/(r, /), it may have another unique positive solution 3j. The condition for 
equation ()22|) to have another solution, that is A (a) 7^ defines the threshold 
above which the infection can establish itself in the population [T%] . 

If we assume that X(0) is proportional to the total population, X(0) = 
bN, where b is a constant, equation (}2*2*j) becomes independent on the popu- 
lation size. In this paper we take b — \i, for simplicity. 
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3.2 Parenterally transmitted HIV 



In this subsection we model the parenteral transmission of HIV due to shar- 
ing contaminated syringes and needles or contaminated blood and blood 
products. 

We consider that sexual transmission among these individuals is negligible 
when compared with the parenterally form of transmission. However, their 
contribution to the sexual spread of the infection will be considered in this 
section. 

Let Yi(a', T)da!dr be the number of individuals parenterally infected with 
HIV (PI), aged between a! and a' + da', who acquired the infection when 
aged between r and r + dr. We consider that Pis enter this group when 
aged d. We also consider that they form a proportion 77 of the population. 
Calling £ the rate of sharing syringes (or receiving blood or blood products) 
multiplied by the probability of getting the infection if the syringe (or blood 
or blood products) is (are) infected, we have: 

y/(a',<r) = X{0)^exp[-f Jl r}exp[-^{T-a l )} 

9(T-a 1 )h 1 (a',T)h 2 (a',r) (23) 

where hi(a', r) and /^(o/, r) were defined in the previous section and are the 
rates of removal from the infective class by death and treatment, respectively. 
Note that we have assumed that in the case of drug users, they remain drug 
users for the rest of their lives and that their mortality rate is not affected by 
the drug addiction. This simplification is partially supported by field work 
we carried out in the past [15] . when we demonstrated that the average time 
of drug usage of the studied community was found to be around 10 years. 
Assuming an antiretroviral treatment rate v (r, a), we have 

Y 2 \\,t,X) = YI{1,t)v{t,\) (24) 

and 

Y 2 I (a,T,l) = Y 2 I (l,r,l)h 3 (a, r, I) 

where /13 (a, r, I) is a function describing the removal of individuals from the 
treated condition by mortality. Then, we get 

Y"/(a, r, /) = hi(a', r)h 2 (a, r)h 3 (a, r, I) u (r, I) (25) 
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Dividing equations and ()25j) by iV, and adding to the corresponding 
terms in equation we finally get, instead of (|T7jl and (J2TJ): 



Kx^r) = X(0)(l-r/)A(r)exp 



X(s)ds — fir 



and 



/ii (a, r) /x 2 (a, r) + r/X(0)£exp [— /xt] 

exp [— £ (r — ai)] 9 (r — a x ) fox (a, r) ft, 2 ( a , r ) 



A(s)<is — /xr 



K 2 (a,r,0 = X(0) (1 - 77) A(r) exp 

hi (I, t) h 2 (I, r) h 3 (a, r, I) v (r, I) 
+i]X(0)^exp [—fir] exp [— £ (r - a x )] 
9 (t - ai) /ii (/, r) /i 2 (Z, t) /i 3 (a, r, Z) z/ (r, /) 
Finally, instead of (J22|) we get 
X(0) 



A(a) 



AT 
exp 



Jo 



(l-^)A(r) 



/ X(s)ds — /jt 
jo 



f(a,a')g{a' - r)hi (a',r) 



A (U) Z 100 /" a 

/x 2 (a', r) drda H V ex P [ — Z 27 "] 

AT Jo Jo 

exp [-£ (r - ai)] 9 (r - a x ) /(a, a')#(a' - r 

X(0) 



/ii (a', r) /i 2 (a', r) drda + 



N 



Jo 



(1-77) 



A(r) exp 



A(s)ds — /xt 



f(a,a')g(a' - r) 



v (r, I) hi (I, t) h 2 (I, t) h 3 (a 1 , r, I) f(a, a') 

X(0) 



A U /"°° /" c 

#i(a' -1,1 - T)dldrda' H — / / 

X Jo Jo 

exp [— /xt] exp [— £ (r — ai)]9 (r — a x ) 



(26) 



(27) 



7/ 



z/ (r, I) hi (I, r) h 2 (I, r) h 3 (a', r, I) 
f(a, a')gi(a' — 1,1 — T)dldrda' (28) 

As mentioned before, if we assume that X(0) is proportional to the pop- 
ulation, for instance, X(0) = fiN, the equation (j2~B*j) becomes independent of 
the population size. 
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4 A schematic model 



In this section we propose forms for the different functions involved in equa- 
tion (|2B|I. 

Let us begin with an untreated population, for which u(t, I) = and 
^2(0', t) = 1 so that the second part of the integral equation (|28|) vanishes. 
As in [3], let us define the rate /(a, a'), in a very schematic form: 

f(a, a') = fo(a, a')9(a — ao)8(a — do) (29) 

where /o( a , a ') is the rate of unprotected sexual contacts and a is the age at 
which individuals enter in the risk behavior group. The Heaviside functions 
8 (a — a ) and 8(a! — a Q ) mean that sexual activity begins after age ao- 

For the function /o(a, a'), which describes the age preferences in acquisi- 
tion of new partners and the age decline in the sexual activity, we propose: 

f (a, a') = Q(3 3 (a)P 3 (a')P,(a - a') (30) 

where 

1 _(>-M) 2 

P 3 (x) = -^e (31) 



Itict 1 

and 

1 



P 4 (a - a') = -=^e CT T (32) 

V 27TCX2 

The parameter M is the age of maximum sexual activity and the constant 
Q is adjusted to give the assumed average number of unprotected sexual 
contacts per unit time. Those forms for the /3-functions (13 1|) and (|32j) were 
chosen for convenience only and, although not supported by social studies, 
they conform with the following facts: a) sexual activity increases with age up 
to a maximum, decreasing thereafter; b) age preferences of an individual are 
distributed around a maximum value which we assumed to occur at the same 
age of the individual. In order to facilitate the calculations we made both 
functions as symmetric around a central value and we assumed a minimum 
age ao, below which there is no sexual activity. In addition we assumed that 
people without treatment, people undergoing treatment, and even people 
with full blown AIDS have the same sexual behavior. This assumption is 
supported by preliminary data in our community of HIV patients (Bueno, 
personal communication) . 
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Let A (a! — r) be the viraemia level of non-treated individuals which 
depends on the time interval since the infection. For g{a! — r) we take the 
form 

g(a' -r) = I(A NT (a' - r)) (33) 

where I(A NT ) is a function representing the transmission of at least one 
infective viral inocula, given the viremic level A NT . As in [3| we assume that 

I {a nt ) = {d + c" log A NT )6 (log A NT - 3) (34) 

where d = —3.35 x and c" = 1.2 xl0 _3 are parameters obtained by 
fitting equation (J34|) to the data by Gray et al. (201 (an almost identical 
relationship between log of viral load and transmissibility was found by using 
data from Vella et al [T7j and Garcia et al [TH], and it is supported by the 
observations of Fideli et al. JHj and Quin et al. jlH])- The Heaviside function 
9 (log A NT — 3) is such that if log A NT < 3 there is no transmission. 

For the removal of infected individuals, we assumed that when (a' — r) < 
L c , where L c is a given critical moment when individuals reach a certain 
viraemia level and are defined as AIDS patients, dies with rate \l. When 
{0! — r) > L c , individuals are subjected to an additional, disease specific, 
death rate, a. A simple form for this removal function is: 



h^a^r) = e-^ a '- T) 6 (L c - (a' - r)) 

+e - M (o'-r)-a((a'-T)-I,c)0 (( a ' _ r ) _ Lc ) (35) 

Let us now consider the treated population. The variables and parameters 
already defined for the untreated population remain the same. We must only 
define values for /i 2 (a / ,r), u(r,l), gi(a' — 1,1 — r) and h 3 (a',r,l), where I is 
the moment in the history of the infection at which individuals begin to be 
treated. 

Let A T (a'—l, l—r) denote the viraemia level after antiretroviral treatment. 
We assumed that the viraemia level after antiretroviral treatment is reduced 
by a certain factor, A(Z — r), increasing thereafter as 

A{a-l,l-r) = (36) 

where ^f(a' — I) is a given function with ^/(0) = 1. Note that with such func- 
tions, the log of the viraemia level just after the beginning of the antiretrovi- 
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ral treatment starts with log A NT (l — r) — log [A(l — r)] and increases with 
log [ty (a' — I)] thereafter. 

Figure 3a schematically illustrates the natural history of HIV infection 
in the absence of antiretroviral treatment. Figure 3b illustrates the effect of 
antiretroviral treatment on the viraemia level. 



Figure 3a 
Figure 3b 



In Figure 3a the viraemia level reaches a maximum immediately after the 
moment r of the infection and drops rapidly after some few weeks, probably 
due to the effect of the immune system. Viraemia then starts to raise, and 
simultaneously a decrease in the counting of CD 4 T-lymphocytes is observed, 
until a critical viral load level is reached (L c ) , when the immune system breaks 
down and full blown AIDS develops. This period lasts for 10 - 15 years. 

Figure 3b represents the assumed natural history of HIV infection in the 
presence of antiretroviral treatment, which begins at the age of infection be- 
tween I and I + dl. The treatment starts at any moment after the viral load 
reaches 30,000 copies/ml [21] (we are well aware, however, that in real life 
other clinical indicators of antiretroviral treatment, like CD 4 counting, are 
used to begin the antiretroviral treatment.). We assumed that the treat- 
ment causes a sharp decrease in the viral load, and that it loses its effect 
immediately so that the viral load starts to rise log-linearly again. This in- 
tends to model the appearance of resistant strains, which we assumed to 
have the same virulence as the original strain. Therefore, the assumption 
of log-linearity in the viral load curve is entirely hypothetical and was in- 
tended to mimic an immediate development of full resistance by HIV to the 
treatment , which is the worst epidemiological scenario. In this situation, we 
are assuming that replication of HIV is no longer constrained by treatment. 
Another assumption of our model is that it takes no account of the role of 
the immunity system on HIV replication. With the above assumptions, the 
viral load after treatment increases until a critical level is eventually reached 
(L' c ), when the immune system breaks down and full blown AIDS develops. 
Hence, the treatment makes the period without AIDS longer than 10 years. 

We assumed also that the antiretroviral treatment schedule is given by 
u(t,1) = v x 9 ((/ — r) — a t ), and it begins, as mentioned above, when vi- 
raemia reaches 30,000 copies/ml, which occurs at the infection age (1—t) = at 
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(~ 6.7 years). Thus 

Since we assume that the treatment affects only the viral load, we may 
define g\(a! — 1,1 — r) as: 

gi (a -1,1-t) = I (A T (a -1,1- r)) (38) 

Finally for h 3 (a,r,l) we may also assume a simple form, similar to 
hiia'-r): 

h 3 (a',r,l) = h 3 (a',l)=e-^ a '^6(L' c -(a'-l)) 

+e -p(a>-l)-8((a'-l)-L> c ) e ^ £/j (gg) 

where L' c is the new critical moment when treated individuals reach a certain 
viraemia level and are defined again as AIDS patients. Those patients are 
subjected to a new differential mortality rate 5 due to AIDS. 

Within this framework we can obtain different particular models by choos- 
ing suitable functions A NT , A T and constants. 

We adopted, for the non-treated patients, the following average values for 
Ai and Lf. 





' Ax (a' 


~ r) 


for 


10 6 

< (a' - r) < Li 




A NT (a' - r) = < 


A 2 (a' 


-r) 


for 


^g3+0.22(a'-r-Li) 

£i < (a' - t) < L c 


(40) 




A 3 (a> 


-r) 




10 6 






- 




for 


(a' - r) > L c 





where L\ = 6 weeks and L c — L\ = years is the period of time it takes 
for the viraemia to reach 10 6 RNA copies per ml since the beginning of the 
second phase. 

Let us now model the effectiveness of antiretroviral treatment by an in- 
stantaneous reduction on the viral load immediately after the introduction of 
the treatment. We consider several effectiveness levels of antiretroviral treat- 
ment, by varying the reduction of the viral load. When the reduction is by 
a factor greater than 10 2 , the antiretroviral treatment is known as HA ART 
(Highly Active AntiRetroviral Treatment). If the reduction is lesser than 10 2 , 
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the antiretroviral treatment is known as non-HAART. Again, for the sake of 
simplicity, we drop the superscript T from the viral load A and introduce a 
subscript j = 4, 5, describing the two possible phases of the natural history 
of the infection after treatment. Note that we assumed no treatment during 
the first phase. 

A T ((a'-l),(l-r)) = 
f A 4 {{a'-l),{l-r)) = lo^K+om-r-L^+o^ia'-i) 

for 0<(a'-l)<L' c 
A 5 ((a'-l),(l-r)) = 10 6 l4ij 
{ for (a'-l)>L' c 

where L' c = 3+k-o.22^-t- Li ) yearg ^ go that _ ^ = 10 o. 22 (o'-0 ; similar to 

the non-treated patients, and A(7 — r) = K > 2 for HAART treatment and 
A(l — r) = K < 2 for non-HAART treatment. In the numerical simulations 
several values of K will be considered. Note that L' c is the period of time it 
takes for the viraemia of treated individuals to reach 10 6 RNA copies per ml 
since the beginning of antiretroviral treatment. 

With the above models of viraemia, g(a' — r) takes the form 



C /(Ai((a'-r))) for < (a' - r) < U 
g{a'-r) = \ J(A 2 ((a'-r))) for L x < (a' - r) < L c , (42) 
[ I(A 3 ((a'-r))) for (a'-r)>L c 

where I(Ai) is a function defined by equation (j2U), representing the proba- 
bility of transmission of at least one infective viral inocula. Correspondingly, 
gi(a' — /, r) takes the form 

I(A 4 ((a'-l),(l-r))) 

for 0<(a'-l)<L' c 
I(A 5 ((a'-l),(l-r))) ^ 
for (a' -l)> L' c 

Finally, for the parenterally transmitted branch of the infection we take 
r] = 0.01, £ = 0.05 per year and a\ = 15 years [T5] , 



9i[a 



U-r) 



5 Numerical results 

In order to analyze the model's performance against an HIV endemic situa- 
tion we solved equation (J28j) numerically for several treatment schedules. For 
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this we divided the population into four classes. Each class obeys an equation 
like equation (|'28|) and intends to mimic a specific risk group, namely, group I 
(GI), with very low level of sexual promiscuity (84% of the total population), 
group II (Gil) with moderate levels of sexual promiscuity (10%), group III 
(GUI) with high levels of sexual promiscuity (5%), and group IV (GIV) with 
very high levels of sexual promiscuity (1%). This somewhat arbitrary divi- 
sion is based on the actual classes of risk recognizable in real populations. 
So GI could represent the general population , Gil promiscuous heterosex- 
uals, GUI male homosexuals, and GIV commercial sex workers and their 
clients.We also assume that in each class there is a fraction rj of Pis, who got 
infected by contaminated syringes/needles or blood/blood products. With 
respect to sexual contacts, it should be stressed that we are only interested 
in unprotected sexual contacts with new partners. In the functions (3 3 and 
(3 4 , which describe the sexual behavior, we set, for all classes, <j\ = 10 years, 
o"2 = 15 years, M = 25 years. The mortality (and fertility) rate was taken 
\i = 1/70 years -1 . The initial age of sexual life a® was set to 15 years. 

The parameter Q is related, in the absence of the infection, to the per 
capita number of unprotected sexual contacts. We can, then, calculate the 
per capita number of unprotected sexual contacts in the population, $, as 

j rco roc 

$ = T7 / / f3(a,a')N{a)da'da 
N Jo Jo 




6{a - a )6(a' - a )N(a')e-> ia da'da (44) 



using equation (J2HJ)- Using equation (fHUJ) we get the relation between the per 
capita number of unprotected sexual contacts in the absence of infection, $, 
and Q: 

/•OO /"OO 

$ = fl 2 Q / (3 3 ( a )(3 3 (a')p 4 (a~a')e-^ a e-^da'da (45) 

Jao Jao 

The parameter Q was adjusted to give the estimated average number of 
unprotected sexual contacts per unit time, $. The results, compatible with 
the literature, for each class were Q = 6 x 10 6 (GI), Q = 1.38 x 10 7 (Gil), 
Q = 1.45 x 10 7 (GUI), Q = 1.55 x 10 7 (GIV), as can be seen in table 1. 

With respect to the PI arm of transmission, as mentioned above, we chose 
the parameters a\ = 15 years, rj = 0.001 and £ = 0.05/year. 

Now, in the presence of HIV infection, without any treatment, the inci- 
dence of sexually transmitted HIV (i) is defined as the number of new cases 
of HIV infection per year per person and was calculated as 
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[i-V) [Jo° A(o) X(q)da] 



AT 

and the incidence of parenterally transmitted HIV (z 7 ) is also defined as the 
number of new cases of HIV infection per year per person and was calculated 
as _ 

i 1 = ?7/i£ / exp [— fj,a] exp [— £(a — ai] 0(a — ai)da (47) 

With the adopted parameters a x = 15 years, 77 = 0.001 and £ = 0.05 per 
year, we have i 1 = 9 x 10 -6 per person-year. 

The prevalence of sexually transmitted HIV (p) in the absence of treat- 
ment is given by 

roo ra 

P= / yi(a,T)drda (48) 
Jo Jo 

where 

yi (a,r) = ^p (49) 

where Yi(a, r) is given by equation (J2fi|) and A(a) given by equation (f2~%J) . with 
z/(r, /) = and ^2(0', t) = 1. The prevalence of parenterally transmitted HIV 
(//), without treatment, is given by 

roc pa 

P 1 = T)ti£exp[-fiT] 
Jo Jo 

exp [— £ (r — ax)] 9 (t — a{) hi(a, r)drda. (50) 

With the adopted parameters a x = 15 years, 77 = 0.001 and £ = 0.05/year, 
we have p 1 = 5.9 x 10~ 4 . 

The results of the simulation for the per capita number of unprotected 
sexual contacts, the prevalence and the incidence of HIV infection at equi- 
librium are given in table 1: 



Table 1: 


prevalence 


and incidence of HIV infect 


ion at equilibrium 






without PI 




with PI 




class(%) 


$ (years -1 


) incidence 


prevalence 


incidence 


prevalence 


GI(84%) 


6.60 |2U 


0.00 


0.00 


8.51E-05 


1.25E-03 


GII(84%) 


15.30 [23 


7.09E-04 


1.08E-02 


1.20E-03 


1.89E-02 


GIII(84%) 


16.00 [2E] 


2.60E-03 


4.50E-02 


2.75E-03 


4.77E-02 


GIV(84%) 


17.10 


4.60E-03 


9.00E-02 


4.63E-03 


9.17E-02 


Average(84%) 


13.20 


2.47E-04 


4.23E-03 


3.75E-04 


6.24E-03 
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The scenario presented in table 1 is not intended to fit real data but rather 
to present qualitatively the effects of treatment. It is interesting, however, 
that the prevalence results are compatible with regions with intermediate 
levels of HIV transmission. By intermediate levels of HIV transmission, we 
mean regions of the world with HIV/AIDS prevalences ranging from 0.5% 
to 1.2% [23]. Note that without the PI class there would be no disease in 
the general population and the infection would be restricted to the other 
three classes. However, as mentioned before, our model is too schematic to 
guarantee that the above conclusion can be extended to any real population. 

Note also that in table 1 the ratio between % and $, which gives the num- 
ber of new cases of HIV infection per unprotected sexual relation with each 
new partner. So, for instance, considering the class GIV, our results point 
to a risk of acquiring HIV infection of approximately 3.0% per unprotected 
sexual contact. This result is of the same order of that reported in the liter- 
ature (for instance, in Thailand this risk was estimated to be between 3.1% 
and 5.6% [23]). 

In order to evaluate the impact of antiretroviral treatment on the inci- 
dence and prevalence of HIV infection we simulated the model for values of 
K (effectiveness of treatment) ranging from to 3 and v (intensity of treat- 
ment) ranging from 0.05 to 0.5. The impact of antiretroviral treatment on 
the incidence and prevalence of HIV infection in each class is shown in figures 
4 to 8. Figure 4 shows the effect of antiretroviral treatment on the class rep- 
resenting GI (Q = 115, 000). In this class the risk of contracting HIV is very 
low, reaching a prevalence without antiretroviral treatment of the order of 
0.1%. In fact, this low prevalence is due to the PI individuals, without whom 
the disease among the general population would disappear, meaning that for 
Q = 115, 000 the infection is below the threshold for its maintenance in the 
general population. Figures 4a and 4b show the incidence and prevalence, 
respectively, as functions of the intensity of antiretroviral treatment v for six 
levels of effectiveness of antiretroviral treatment K. Note that the infection 
drops monotonically with v for all values of K. 

Figure 4a 
Figure 4b 

Figure 5 shows the effect of antiretroviral treatment on the class rep- 
resenting the subpopulation Gil (Q = 138,000). In this class the risk of 
contracting HIV is higher than in the GI, described above. Without an- 
tiretroviral treatment the prevalence of HIV infection reaches the order of 
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2.0%. Figures 5a and 5b show the incidence and prevalence, respectively, as 
functions of the intensity of antiretroviral treatment v for six levels of effec- 
tiveness of antiretroviral treatment K. Note that for K = 0.5 the treatment 
results in higher prevalence and incidence than in its absence, for all values of 
v. Moreover, both the incidence and the prevalence show an initial increase 
with the treatment intensity and a decrease after around v = 0.15. 

Figure 5a 
Figure 5b 

To understand the phenomenon described above, one should consider the 
following: first, a low intensity treatment means that the individuals, on 
the average, start to be treated later than with a high intensity treatment; 
secondly, starting the treatment implies in reducing the probability of HIV 
transmission proportionally to the reduction in the log of viral load due to the 
effect of the treatment; finally, treating individuals implies in a longer survival 
period, therefore increasing the total number of sexual contacts of those 
individuals. Hence, if the effectiveness of the treatment (K) does not reduce 
the log of the viral load sufficiently to decrease the probability of transmission 
per sexual contact such as to compensate the increased transmission of HIV 
due to the higher number of sexual contacts, the total contribution of those 
individuals to HIV transmission will increase. Otherwise it will decrease. 
This effect repeats itself for the other simulations described below. 

Figure 6 shows the effect of antiretroviral treatment on the class repre- 
senting the subpopulation of GUI individuals (Q = 145,000). In this class 
the risk of contracting HIV is higher than that of Gil, described above. 
Without antiretroviral treatment the prevalence of HIV infection reaches the 
order of 4.8%. Figures 6a and 6b show the incidence and prevalence, respec- 
tively, as functions of the intensity of antiretroviral treatment v for six levels 
of effectiveness of antiretroviral treatment K. Note that, for K = 0.5 the 
prevalence and the incidence is greater than without treatment for all values 
of v. For K > 1, the incidence drops monotonically with v. However, for 
K — 1, the prevalence is greater than that without treatment for values of v 
up to approximately 0.8. 

Figure 6a 
Figure 6b 
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Figure 7 shows the effect of antiretroviral treatment on the class rep- 
resenting the subpopulation GIV (Q = 155,000). In this class the risk of 
contracting HIV is the highest among all classes considered. Without an- 
tiretroviral treatment the prevalence of HIV infection reaches the order of 
9%. Figures 7a and 7b show the incidence and prevalence, respectively, as 
functions of the intensity of antiretroviral treatment v for six levels of ef- 
fectiveness of antiretroviral treatment K. Note that both the incidence and 
the prevalence for K = 0.5 are higher than those without treatment, for all 
values of v. For K — 1 and K = 1.5, the incidence drops monotonically but 
the prevalence is greater than that without treatment for small values of v. 

Figure 7a 
Figure 7b 

Figure 8 shows the weighted average of the incidence (8a) and prevalence 
(8b) curves over the entire population. Note that both the incidence and the 
prevalence for K = 0.5 are higher than those without treatment, for all values 
of v. For K > 1, both the incidence and prevalence drop monotonically. 

Figure 8a 
Figure 8b 

6 Comments and conclusions 

In this paper we presented a very simple model of the steady-state effect of 
HAART on HIV incidence and prevalence. The model mimics current treat- 
ment guidelines applied in Brazil. However, the model does not intend to fit 
the data with any acceptable degree of accuracy since detailed information 
necessary are not available. So, this paper intends to provide a conceptual 
and mechanistic understanding of the possible long term effects of treatment 
on the dynamics of HIV transmission. As mentioned by Anderson }29| . one 
of the purposes of modelling is to help identifying areas in which better epi- 
demiological data is required to refine prediction and improve understanding, 
guiding, in a way, field research. We hope our model can be useful in pointing 
which parameters should be better determined in future studies. 

The model could be extended to allow the calculation of temporal evo- 
lution of the effects of treatment on HIV incidence, as we did for rubella 
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vaccination in [T2|. At the present stage we calculated only steady-state ef- 
fects, which is important enough for assessing long term effects of treatment 
on trends of HIV dynamics. 

One important aspect of the model is that the population was divided 
into four compartments according to their sexual activities. We also assumed 
another compartment of Pis and we assumed that only they interact with the 
other four different model compartments. The Pis were singled out because 
they acquire the infection by a different route. It would appear natural to 
consider that the risk groups would interact with each other. However, this 
would introduce more unknown parameters and, therefore, for simplicity we 
considered such an interaction as negligible. Note that this implies that if 
an individual from any single risk group, for instance GI, has an unprotected 
sexual contact with any other individuals from another group he/she and 
his/her stable partner would be considered as belonging to this class. 

Our main conclusions are as follows: 

1. The disease, according to the model, is almost completely wiped out 
when we consider the most effective treatment (K = 3.0) simulated. 
This conclusion should be taken with great care since it may be the re- 
sult of some 'mathematical pathology' taken to its extreme. In fact, the 
model predicts that under this treatment regime the disease is main- 
tained in the population (all four classes) due to the interaction of indi- 
viduals from the other classes with Pis. Again, this conclusion should 
be taken with caution since, by allowing strong interactions between the 
distinct sexual behavior classes the disease might not disappear under 
treatment. In addition, note that our simulated treatment represents 
the worst scenario (the immediate evolution of complete resistance with 
no further alteration in the treatment scheme), while in clinical practice 
modifications of treatment schemes should always follows a significant 
increase in viraemia. 

2. The impact of the treatment on HIV incidence or prevalence depends 
on the level of sexual activity of the subpopulations considered, being 
more pronounced on the subpopulations with the highest sexual ac- 
tivity levels. By impact of treatment we mean the difference between 
the pre-treatment level of incidence or prevalence and the equilibrium 
attained with the maximum intensity of treatment, v. This conclusion 
is valid only for the most effective treatment scheme, K = 3.0 
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3. Inefficient treatment, K < 1.0, can be prejudicial on the subpopulations 
with high levels of sexual activity. For instance, in populations with 
intermediate levels of sexual activity ($ = 0.153 years -1 ), the effect 
of inefficient treatment depends on the intensity of treatment, z/, in a 
curious way. For v between and 0.15, there is an increase in both 
incidence and prevalence, which then drops thereafter. The reason for 
this behavior is explained in the main body of the text. 

As a general comment, there are many ways to express the intensity of 
transmission of an infection, among which the classical basic reproduction 
ratio, Rq, the force of infection, A(a), the incidence and the prevalence. The 
basic reproduction ratio is the greatest eigenvalue of the Frechet derivative 
with respect to A (a) of the operator, which is the right hand side of equation 
I2HJ calculated at A (a) = (see [12], [SOD- In this paper we calculated the 
incidence and the prevalence of HIV as these are the parameters most used 
by public health authorities to monitor HIV epidemic. 

Finally, since HIV treatment begins late in the second phase (in practice 
and in this model), the intensity and diversity of effects that result from 
different treatment strategies is another evidence of the importance of the 
asymptomatic phase of HIV infection on the spread of the virus. This has 
already been pointed out pi , contrasting with former opinions of some in- 
fectious disease practitioners and epidemiologists. 
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Captions for the figures 



Figure 1: Reduction in the mortality by AIDS in the period between 
1996 and 2001. The total number of averted deaths summed up to 90,000 
patients in this period. 

Figure 2: Reduction in the incidence rates in Brazil two years after the 
introduction of universal highly active antiretroviral treatment. 

Figure 3: (a) Model assumed to describe the natural variation of HIV 
viraemia along the natural history of the infection (shaded area). The infec- 
tion happens at age between r and r +dr and L c marks the beginning of full 
blown AIDS, (b) Model assumed to describe the natural variation of HIV 
viraemia along the natural history of the infection in presence of antiretro- 
viral treatment (shaded area). The infection happens at age between r and 
r +dr, the treatment begins at age between/ and I + dl and L' c marks the 
beginning of full blown AIDS. 

Figure 4: (a) The incidence of HIV in the GI (general population) group 
as a function of the treatment intensity rate v for several levels of effectiveness 
of antiretroviral treatment K. (b) The prevalence of HIV in the GI (general 
population) group as a function of the treatment intensity rate v for several 
levels of effectiveness of antiretroviral treatment K. 

Figure 5: (a) The incidence of HIV in the class representing the sub- 
population of Gil (promiscuous heterosexuals) group, as a function of the 
treatment intensity rate v for several levels of effectiveness of antiretroviral 
treatment K. (b) The prevalence of HIV in the class representing the sub- 
population Gil (promiscuous heterosexuals), as a function of the treatment 
intensity rate v for several levels of effectiveness of antiretroviral treatment 
K. 

Figure 6: (a) The incidence of HIV in the class representing the sub- 
population GUI (male homosexuals), as a function of the treatment intensity 
rate v for several levels of effectiveness of antiretroviral treatment K. (b) The 
prevalence of HIV in the class representing the sub-population GUI (male 
homosexuals), as a function of the treatment intensity rate v for several levels 
of effectiveness of antiretroviral treatment K. 

Figure 7: (a) The incidence of HIV in the class representing the sub- 
population GIV (sex workers) and their clients, as a function of the treatment 
intensity rate v for several levels of effectiveness of antiretroviral treatment 
K. (b) The prevalence of HIV in the class representing the sub-population 
GIV (sex workers) and their clients, as a function of the treatment intensity 
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rate v for several levels of effectiveness of antiretroviral treatment K. 

Figure 8: (a) The incidence of HIV in the population average as a func- 
tion of the treatment intensity rate v for several levels of effectiveness of 
antiretroviral treatment K. (b) The prevalence of HIV in the population 
average as a function of the treatment intensity rate v for several levels of 
effectiveness of antiretroviral treatment K. 
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Figure 1 
AIDS mortality in Brazil 




Figure 2 
Incidence of AIDS in Brazil 
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